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ABSTRACT 

Be stars possess gaseous circumstellar decretion disks, which are well described using standard a-disk the- 
ory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during 
which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous de- 
cretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter 
a = 1 .0 ± 0.2, corresponding to a mass injection rate M = (3.5 ± 1 .3) x 10~ 8 M yr" 1 . Such a large value of a 
suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock 
dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass loss rate 
inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, 
but some other mechanism. 

Subject headings: circumstellar matter — stars: individual (28 CMa) — stars: activity — stars: emission-line, 
Be 



1. INTRODUCTION 

It is now generally accepted that Classical Be stars pos- 
sess gaseous circumstellar disks that ar e responsible for 
their excess emission (fo r reviews, see iPorter & Riviniusl 
12003b iBaade et all 1201 ll) . Furthermore spectroastrometry 
and optical interferometic measure ments indicate that this 
material is in K eplerian rotation (lOudmaiier et al.l 1201 ll 
iRraus et al. 2011). Consequently, as detailed in the review 
of ICarcioril (1201 If), th e steady-state viscous decretion disk 
model dLee et alii 199 lb is capable of explaining many of the 
(time-averaged) observational properties of Be star disks. 

Most likely these disks are fed by mass-loss from 

their rapidly rotating central star s. Both observational 

jRivinius etafl Il998t iNeiner et al.l |2002|) and theoretical 
( Ando 119861: ICranmerl 120091) studies suggest that non-radial 
pulsations may play an important role. However, to create a 
viscous decretion disk, this material must be placed into or- 
bit. Since the stars most likely are rotating at less than their 
critical speed (e.g.. lCranm er 2005), the mass loss mechanism 
must add some angular momentum before injecting the ma- 
terial into to the disk. From a theoretical point of view, this 
requirement is problematic, so, as ye t, no entirely s atisfactory 
mechanism has been proposed (see Owocki 2006], for a de- 
tailed assessment). 

Regardless of how material is injected into the disk, if this 
material is in orbit, turbulent viscosity will redistribute it on 
a viscous diffusion timescale, producing an outflowing decre- 
tion disk, as long as the disk is fed continually. However, Be 
stars exhibit a wide range of variability, ranging from emis- 
sion line profile variations (e.g, V/R variability, indicative of 
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densi ty waves within the disk; Ste fl et al.l2009l:ICarciofi etakl 

2009) , to photometric and polarimetric variability, indicating 
that frequently the disks are not fed at a constant rate. In some 
cases, the disk c an disappear for decades b e fore being rebuilt 
(see r eview by lUnderhill & Doazanlll982t IWisniews ki et al.l 

2010) . Oftentimes, this rebuilding phase occurs in a series of 
outbursts, during which the disk apparently is fed for a short 
period of time followed by a period of quiescence, during 
which the nascent disk can partially reaccrete onto the star. 
Similar outbursts can occur while the disk is disappearing. 

Once such example is the Be star 28 CMa. Over the past 
40 years, 28 CMa has exhibited quasi-regular outbursts, ev- 
ery 8 years or so, when the star brightens by about half a 
magnitude in the V-band. This usually occurs via multiple 
individual injections (e.g., the rapid rises seen in Fig. 1), dur- 
ing which the star brightens at a rate of 0.02±0.01magd _1 , 
due to free-bound and free-free radiatio n from the disk. The 
outbursts last for about two years (see Ste fl et al.l 120031 for 
a detailed discussion of the long-term photometric variations 
of 28 CMa). A particularly well-monitored outburst occurred 
from 2001 to 2004, and this outburst was followed by a long 
quiescent phase from 2004 to 2009. Furthermore it appears 
that during this quiescent phase, very little additional material 
was injected into the disk (only a few minor events). Thus na- 
ture has provided us the perfect experiment to study how Be 
star disks dissipate. 

If we accept as a premise that the viscous decretion disk 
model is correct (and that no other mechanisms participate in 
the disk mass budget), then the timescale of the post-outburst 
photometric decline is controlled solely by the disk viscosity. 
In this paper, we model the V-band light curve with the goals 
of: 1) directly measuring the disk viscosity parameter, a, and 
2) further testing the viscous decretion disk model by com- 
paring its time-dependent predictions of the disk dissipation 
against the detailed shape of the photometric light curve. 

2. OBSERVATIONS 

Visual photometric observations of 28 CMa (cj CMa, 
HD 56139, HR2749; B2-3IV-Ve; m v = 3.6 to 4.2, B-V = 
-0.1 to -0.2) were made from early 2003 t o late 2009 using 
a modified version of the Argelandej] d 1 8431) method in which 
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FIG . 1 . — V-band Light Curve. Visual observations of 28 CMa (grey triangles) are shown in comparison to our model fits for different values of a, as indicated. 
Phase I (MJD = 52900^53070) is the initial decline, which was used to determine the value of a. Phase II (MJD = 53070-54670) is the slow disk-draining phase. 
The inset shows the reduced chi-squared of the Phase I fit for different values of a. The horizontal dotted line indicates the 90% confidence level. 

diative transfer calculation of the emergent flux at selected 
times. The disk s urface density is mode led with a vis- 
cous decretion disk dLee. Osaki. & S aio 1991) using the time- 
dependent hydrodyn amics code singlebe (lOkazaki et alj 
2002; Okazaki 2007), which permits time-variable mass in- 
jection at a point placed just outside the stellar surface. This 
code solves the isothermal 1-D time-dependent fluid equa- 
tions (iLvnden-Bell & Prin^ldfl 974t iPringldl 1981b in the thin 
disk approximation. The output is the disk surface density, 
£(r, t ), as a function of radius, r, at selected times, t = t„ as 
well as the mass flow rate as a function of position M(r). This 
surface density is converted to volume density using the usual 
vertical hydrostatic equilibrium solution (a Gaussian) with a 
power law scale height, H = ho(r/R) 15 , where R is the stel- 
lar equatorial radius. The volume density then is used as the 
input to our 3-D Monte Carl o radiative transfer code, HDUST 
(ICarciofi & Bjorkmanl 120061) . HDUST simultaneously solves 
the NLTE statistical equilibrium rate equations and radiative 
equilibrium equation to obtain the hydrogen level popula- 
tions, ionization fraction, and electron temperature as a func- 
tion of position. Its output is the emergent SED (and other 
observables of interest), which is used to calculate the V-band 
excess, AV, of the disk as a function of time. We assume that 
this variable disk emission dominates any intrinsic variability 
of the star. 

We determine the V-band excess of 28 CMa at the on- 
set of the steep decline in late 2003 (time ?o) by calculat- 
ing the average pre-decline brightness (my = 3.63 ±0.05) and 
subtracting the average near the end of the quiescent phase 
{my = 4.18 ±0.02). This gives an initial excess AV(fo) = 
-0.54 ± 0.05, where the error is estimated from the standard 
deviation of the means. Although the disk building typically 
is not steady, the growth time required to produce this entire 
excess f g = 0.54 mag/0.02 mag d -1 = 0.07 yr. 

Since the brightness of 28 CMa was relatively steady for 
two years prior to the decline, we assume the initial sur- 
face density of the di sk was nearly the steady-sta te solution, 
£(r,f ) = Zo(r/Ry 2 dBjorkman & Carcioffe OOS). To match 



TABLE 1 
Stellar Parameters 

Parameter Value 

L 5224L© 

r Hi , 22000 K 

Rpole 6.0 Rq 

loggpole.-. 3.84 

M 9 M Q 

V m , 350 km s 



-i 



y cril 436kms-' 

Req 7.5 Rq 

i 15° 

the estimation of the variable star's brightness is performed 
by comparison with two non-variable stars, one brighter, the 
other fainter than the target. The accuracy depends on the 
availability of nearby comparison stars with similar bright- 
ness and color. Since 28 CMa is variable, several pairs of 
comparison stars are required. Based on comparisons with 
simultaneous photoelectric V- band photometry, the error is 
typically less than .05 mag (lOtero. Fraser. & Lloydl 1200 ll : 
lOtero & Moonll2006l) . However, based on the typical disper- 
sion of the 28 CMa observations, we adopt a uniform 1-er error 
of 0.03 mag. The resulting V-band light curve for 28 CMa is 
shown in Figure Q] 

3. MODELING 

To model the circumstellar disk of 28 CMa, we require 
physical parameters for the central star, including th e incli- 
nation angle of the system. Fortunately. iMaintz et al.l (2003) 
performed an extensive investigation of the non-radial pulsa- 
tion modes of this star. Using detailed fits of the line profile 
variability, they conclude that the star is viewed nearly pole- 
on at an inclination angle i = 15°. Table Q] lists their recom- 
mended values for the stellar parameters. 

Given the stellar parameters, we now model the disk dissi- 
pation of 28 CMa. To do so, we combine a time-dependent 
calculation of the disk surface density with a NLTE ra- 
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FIG. 2. — Disk dissipation. Shown is the V-band excess as a function of 
time after mass injection stops for different values of a as indicated. 
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0.8 
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2.9 


2.9 
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the initial AV, we require So = 1.6 ± 0.5 gem -2 , whose er- 
ror is dominated by the uncertainty in the inclination angle 
(i = 0°-30°). This implies a pre-decline disk mass Ma = 
2irJ: R 2 ln(r/R) = 3 x 10" 9 M Q interior to r = 10R. 

To produce the initial steady-state disk, we run SINGLEBE 
for a long time with a constant mass injection rate, M. Note 
that Eo depends on the ratio M/a, so whenever we change the 
disk viscosity parameter, a, we adjust M to maintain the de- 
sired value for £rj. After SINGLEBE achieves steady-state, we 
turn off the mass injection (by setting M = 0) and let the disk 
dissipate. Figure [2] shows the decline in the V-band excess 
as a function of time. Note that as the viscosity parameter 
a increases, the disk dissipates more rapidly. Initially there 
is a rapid decline as the inner disk (where the V-band excess 
is produced) adjusts from outflow to infall and is reaccreted. 
Subsequently, the remaining mass of the disk is drained via 
quasi-steady-state accretion, and the V-band excess continues 
to slowly drop to zero. 

Given the model curves in Figure [2] we fit the Phase I por- 
tion of the observed light curve (see Fig. [TJ by adjusting only 
two parameters: a and to- Table [2] lists the best-fitting value 
of to for each value of a and the reduced chi-squared value 
for the Phase I data, x\ ■ Additionally, we list the reduced chi- 
squared values of the extrapolation of the fit through Phase II, 
xji, as well as the reduced chi-squared for the entire data set, 
X 2 - Finally, each model curve is shown in comparison to the 
observations in Figure Q] The overall best fit to the Phase I 
data is a = 1.0 ±0.2 (90% confidence interval). Note that 
when we extrapolate the fit to Phase II, the fit is still reason- 
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FIG. 3. — Disk Structure. Shown are the disk temperature (top panel), 
hydrogen ground state population (middle panel), and density (bottom panel) 
as a function of radius for the initial state of the disk (solid line), at the end 
of the Phase I (dotted line) and at the end of the simulation (dashed line) for 
the a = 1 .0 case. The circle marks the position of the stagnation point (the 
division between inflow and outflow) for the 2004. 1 curve. 

able, requiring a large value of a, albeit with a larger uncer- 
tainly. 

4. DISCUSSION 

As may be seen in Figure Q] our model fits quite well the 
detailed shape of the V-band decline observed in 28 CMa. 
Since the V-band excess is pr oduced in the i nnermost parts 
of the disk (typically r < 2R, ICarciofi 1 120 1 lh . we infer that 
the viscous decretion disk model correctly predicts the time- 
dependent evolution of the disk density (at least at small 
radii). Not only does the model reproduce the rapid early 
decline (Phase I), it also reproduces the slow disk-draining 
(Phase II) with a single physical model. Furthermore, there 
is no evidence that the value of a changed from Phase I to 
Phase II, despite the rather large change in physical condi- 
tions in the disk shown in Figure [3] 

When the mass injection ceases, the disk readjusts most 
quickly at small radii. This is because the viscous diffusion 
timescale, t^g oc r 1//2 a _I , increases with radius. Reaccretion 
of material from the inner disk occurs because the turbulent 
viscosity transports angular momentum outward, transferring 
angular momentum from the inner disk to the outer. Hence, 
the outer parts of the disk continue to expand while the inner 
parts move inward. This produces simultaneous inflow in the 
inner disk and outflow in the outer disk. As a function of time, 
the stagnation point (the division between inflow and outflow) 
moves outward, controlled by the diffusion time. Interior to 
the stagnation point, the disk adopts a quasi-steady accretion 
solution, while exterior to the stagnation point the disk re- 
mains in quasi-steady decretion (see Fig. [3). Hence the outer 
disk acts as a mass reservoir that feeds reaccreting material 
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through the inner disk where the V-band excess is produced. 

We conclude that the decline rate in Phase I is controlled by 
the viscous diffusion time; this is what determines the value of 
a in our model fit. Since a is the only free physical parameter 
that we can vary, the fact that we can independently reproduce 
the rate of decline in Phase II (when material is being trans- 
ferred from the outer through the inner disk) is a non-trivial 
test of the time-dependent viscous decretion disk model. 

In our analysis we assume, during the dissipation, that the 
disk is not being fed and no other important mass loss oc- 
curs. Two such mass loss mechanisms are disk ablation and 
entrainment by the stellar wind, but general estimates o f their 
mass loss rates dKrticka et alJl201 ltlGayley et al.ll 19991) indi- 
cate they are much smaller than that of the stellar wind. The 
disk reaccretion rate during Phase I is much larger than the 
wind mass loss rate, owing to the large value of a, conse- 
quently neither plays an important role. 

Perhaps the most surprising aspect of our results is the ex- 
tremely large viscosit y parameter (a sa 1 .0) tha t is required. 
In the a-disk model (Shakura & Sunyae~vl ll973l) . a parame- 
terizes the turbulent (eddy) viscosity, v = v tur b^ = ac s H, where 
Vturb is the flow speed of the eddies and I is their size scale. 
Generally, one assumes that supersonic turbulence will not oc- 
cur because of shock dissipation, so Shakura & Sunyaev sup- 
pose the maximum flow speed is the sound speed, c s , and that 
the characteristic size of the circulation cells is the disk pres- 
sure scale height, H. Hence one expects < a < 1 . To pro- 
duce a = 1 .0 requires either roughly sonic flow and/or some- 
what larger circulation cells. If in fact v tul -b « c s , this would 
strongly suggest that the origin of the turbulent mixing is an 
instability within the disk that grows until shock dissipation 
halts the growth. Such an instability would also explain why 
a « 1 regardless of the physical conditions within the disk. 

Having now directly measured the value of a from the 
V-band decline rate, we can determine the mass decretion 
rate of 28 CMa. Normally one cannot use the disk den- 
sity to determine the mass loss rate because the radial flow 
speeds are too s mall to be measured from the line profiles 
dHanu schik 2000). Similarly, from a theoretical point of view, 
the disk surface density depends on the ratio M/a, so one 
must know a to find M. Assuming the pre-decline value 



of a is the same as that we measure during disk dissipation 
(a = 1 .0), we find that the mass injection rate of 28 CMa 
was M = (3.5 ± 1.3) x 10" 8 M Q yr" 1 . An important check 
of this rate can be made using the growth time, f g , in com- 
bination with the pre-decline disk mass, Mj, which implies 
M ~ M d /fg = 4 x 10" 8 M Q yr" 1 during the growth phase. 
Consequently, we infer there is no evidence for a significant 
change in a between the disk growth and the disk dissipation 
phases. 

We find, therefore, that the mass injection rate is much 
larger than previously thought, and is at least an order of mag- 
nitude larger than the ma ss loss in the stellar wind, as inferred 
from UV observations dSnowl[l981l) . This implies that the 
stellar wind is not the mechanism responsible for mass in- 
jection into the disk, unless the wind can be prevented from 
producing strong, easily observable UV wind lines. 

5. CONCLUSIONS 

Light curve fitting to the dissipation of a disk can be a pow- 
erful tool for measuring a. Further it provides additional 
quantitative tests of the viscous decretion disk model. Here 
we have shown that the viscous decretion disk model can re- 
produce the detailed time-dependence of the dissipation of a 
Be star disk. Using the V-band excess to measure how the 
disk density decreases with time, we have been able for the 
first time to directly measure the value of the disk viscosity 
parameter. We find a = 1.0 ±0.2, which provides an impor- 
tant clue about the origin of the turbulent viscosity, suggesting 
that it likely is produced by an instability in the disk whose 
growth is limited by shock dissipation. 

Finally, knowing the value of a, we calculate that the mass 
decretion rate of 28 CMa was M = (3 .5 ± 1 .3) x 10" 8 M© yr" 1 . 
Such a large required mass injection rate places strong con- 
straints on any proposed stellar mass loss mechanism. 



We thank the referee, Stanley P. Owocki, for his use- 
ful comments. We acknowledge support from CNPq grant 
308985/2009-5 (ACC) and Fapesp grants 2010/19029-0 
(ACC), 2010/16037-2 (JEB), and 2009/07477-1 (XH). 



REFERENCES 



Ando. H. 1986, A&A, 163, 97 

Argelander, EW.A. 1843, Uranometria Nova, Berlin 

Baade, D., Rivinius, T., Stefl, S., & Martayan, C. 2011, IAU Symposium, 
272, 1 

Bjorkman, JE. & Carciofi, A.C. 2005, in ASP Conf Ser. 337, The Nature 
and Evolution of Disks Around Hot Stars, ed R. Ignace & K. Gayley (San 
Francisco: ASP), 75 

Cai'ciofi, A. C. & Bjorkman, J. E. 2006, ApJ, 639, 1081 

Carciofi, A. C, Okazaki, A. T., Le Bouquin, J.-B., et al. 2009, A&A, 504, 
915 

Carciofi, A. C. 201 1, IAU Symposium, 272, 325 
Cranmer, S. R. 2005, ApJ, 634, 585 
Cranmer, S. R. 2009, ApJ, 701, 396 

Doazan, V., Franco, M., Rusconi, L., Sedmak, G., & Stalio, R. 1983, A&A, 
128, 171 

Gayley, K. G., Owocki, S. P., & Cranmer, S. R. 1999, ApJ, 513, 442 
Hanuschik, R. W. 2000, IAU Colloq. 175: The Be Phenomenon in Early-Type 

Stars, 214, 518 

Kraus, S., Monnier, J. D., Che, X., et al. 201 l, larXiv:110 9.3447 

Krticka, J., Owocki, S. P., & Meynet, G. 2011, A&A, 527, A84 

Lee, U., Osaki, Y, & Saio, H. 1991, MNRAS, 250, 432 

Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 

Maintz, M. Rivinius, Th., Stefl, S., Baade, D., Wolf, B., & Townsend, R. H. D. 

2003, A&A, 411, 181 
Neiner, C. 2002, A&A, 388, 899 



Oudmaijer, R. D., Wheelwright, H. E., Carciofi, A. C, Bjorkman, J. E., & 

Bjorkman, K. S. 201 1, IAU Symposium, 272, 418 
Okazaki, A. T, Bate, M. R., Ogilvie, G. I., & Pringle, J. E. 2002, MNRAS, 

337, 967 

Okazaki, A. T. 2007, Active OB-Stars: Laboratories for Stellar and 

Circumstellar Physics, 361, 230 
Otero, S. A., Fraser, B., & Lloyd, C. 2001, Information Bulletin on Variable 

Stars, 5026, 1 

Otero, S. A., & Moon, T. 2006, The Journal of the American Association of 

Variable Star Observers, 34, 156 
Owocki, S. 2006, Stars with the B[e] Phenomenon, 355, 219 
Porter, J.M. & Rivinius, T. 2003, PASP, 115, 1153 
Pringle, J. E. 1981, ARA&A, 19, 137 

Rivinius, Th., Baade, D., Stefl, S., Stahl, O, Wolf, B., & Kaufer, A. 1998, 

A&A, 333, 125 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Snow, T. P., Jr. 1981, ApJ, 251, 139 

Stefl, S., Baade, D., Rivinius, T, et al. 2003, A&A, 402, 253 
Stefl, S., Rivinius, T., Carciofi, A. C, et al. 2009, A&A, 504, 929 
Underhill, A. & Doazan, V. eds. 1982, B Stars with and without Emission 

Lines (NASA SP-456; Washington, DC: NASA) 
Wisniewski, J. P., Draper, Z. H., Bjorkman, K. S., et al. 2010, ApJ, 709, 1306 



